Auxetic properties of polycrystals 
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Young's and shear moduli and Poisson's ratio of polycrystaUine solids consisting of 2D quadratic 
and 3D cubic randomly oriented grains of the same size and shape is studied. Considered polycrystals 
are initially unstrained. It is shown that for such polycrystals the division of the mechanical stability 
regions into areas of various auxeticity properties is different than for monocrystalline solids. In 
particular the regions of complete auxeticity enlarge. 
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I. INTRODUCTION 

In our recent papers we considered mechanical char- 
acteristics of all 2D symmetry systems [l| and of 3D 
crystalline structures of high and middle symmetry 
In particular we derived explicit expressions for Young's 
E and shear G moduli as well as for Poisson's ratio v 
depending on directional cosines of angles between di- 
rections n of the load and the direction m of the lateral 
strain, respectively, and directions of the crystalline sym- 
metry axes. 

Our particular attention was paid to 2D and 3D 
cubic [J] structures. All initially unstrained quadratic 
and cubic materials are characterized by three parame- 
ters si, S2 and S3 belonging to a half-infinite [s\ > 0) 
prisms with (stability) triangles (ST) in the base lying in 
the (s2,S3)-plane. 

Among mechanical characteristics, Poisson's ratio v 
is particularly interesting because the phenomenon of a 
negative i.e. solids expanding transversely to an ap- 
plied tensile stress, is counterintuitive. Elastic materials 
with negative v are termed auxetics. Materials with neg- 
ative v for all pairs of vectors n, m are complete aux- 
etics. Materials are auxetics if there exists pairs of n, 
m for which v < 0. Materials are non-auxetics if v > 
for all pairs of vectors n, m. In papers [|[ and [3] we 
established regions of ST in which crystalline quadratic 
and cubic materials are completely auxetic, auxetic and 
non-auxetic. One may expect that complete auxetics are 
particulary promising from the standpoint of technology. 

Many technically important materials are polycrys- 
taUine. In view of potential applications of crystalline 
auxetic materials, the question how polycrystallinity in- 
fluences the auxetic properties deserves some attention. 

The elastic properties of polycrystaUine materials de- 
pend on the single-crystal elastic constants of the crystal- 
lites which build up the polycrystal and on the manner in 
which crystallites are connected. In most cases exact ori- 
entations, shapes and connections are not known. These 
gaps in knowledge can be overcame with the help of orien- 



tation and grain shape distribution functions. Generally, 
one has to resort to more or less realistic assumptions 
(cf. (H). We assume that the grains have the same shape 
and volume and that their orientations are completely 
random. 



II. MECHANICAL CHARACTERISTICS OF 
MONOCRYSTALLINE 3D AND 2D CUBIC 
ELASTIC MATERIALS 

A. Relations between the stiffnesses and 
compliances and components of vectors n and m 

Mechanical properties of crystalline materials depend 
on the components of the stiffness tensor C or on the 
components of compliance tensor S. In the coordinate 
system related to the crystalline axes these tensors have 
the minimal number of independent components that de- 
pend on the particular crystalline symmetry. In the case 
of cubic and quadratic (i.e. 2D cubic) materials, these in- 
dependent components are Cn, C12 and Cee and Su, S'12 
and See respectively (cf. Q). The components C/j 
and Su (I, J = 1, 2, 3, . . . , 6) are related. For quadratic 
materials [l[ 
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The direction n of applied tension and direction m in 
which the lateral expansion/contraction is measured are 
mutually perpendicular, i.e. nm = 0. For quadratic 
materials 
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n 2 : 



m 2 = cos ip, 
—mi = sin ip. 



(3) 



2 



In the case of cubic materials [9j 

n\ = cos a cos tp cos — sin a sin 9, 
n 2 = cos a cos <p sin + sin a cos 0, 
ri3 = — cos a sin <p, 



(4) 



and 



where 



mi = smt/Jcos w, 
m 2 = sin^sin^, 
m3 = cos ip, 

TT<a<TT,0<ip<TT,0<8< 2n. 



(5) 
(6) 



B. Mechanical characteristics of 2D and 3D cubic 
materials 

In the case of quadratic materials we obtained [l[ 
E-\n) = sj/2 + s L (nl-nD/2 



+2sMn\n\ > 0, 



(7) 



[4G g (m, n)] 1 = s L (mini - m 2 n 2 ) 2 /2 



^(m, n) 
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+% (miti 2 + m 2 ni) 2 /2 > 0, (8) 

s.//2 + s L (m? - mj) - n^) /2 
+2s M TOiTO 2 n 1 n 2 . 

(9) 

In the case of cubic materials in paper we obtained 
the familiar formulas [|[ 

E- 1 (n) = [(sj-s M )/3 + s L } 

+ (s M -s L )p(n), (10) 
^c(m, n) 

~^(nT = (Sj M)/3 

+ (s M -s L )P(m,n), (11) 

[4G c (m, n)]- 1 = s L /2 + (s M - s L ) P(m, n), (12) 
where 

P(n)=^nf, (13) 



P (m, n) = P (n, m) = ^ (m,ti,)' 



(14) 



The coefficients sj, sl and sjj are the eigenvalues 
of the compliance tensor S. They are related to the 
eigenvalues of the stiffness tensor C, namely s u = c^ 1 
(u — J, L, M). For quadratic materials [l[ 

s.j = Sn + S12, s L = (Su - Si 2 ) , s M = 5*66/2. (15) 

For cubic materials Q 

s.j = Sn + 2S12, s L = S 66 /2, s M = Sn - S 12 . (16) 



III. THE STABILITY CONDITIONS FOR 
CUBIC AND QUADRATIC MATERIALS 

Both stiffness and compliance tensors are positive, 
therefore their eigenvalues su, cjj (U — J, L, M) are also 
positive. This yields the familiar inequalities (cf. @, [l|), 
namely 

An, A m > 0, {An - A 12 ) > (A = G, S), (17) 

and 

A u + {d-l)A 12 >0, (18) 

where d depends on dimensionality D: d = 2 for 2D 
and d = 3 for 3D. We shall underline that the above 
inequalities hold only for initially unstrained crystals [f| . 

For quadratic materials more reasonable choice of in- 
dependent parameters was proposed in our paper , and 
for cubic materials by Every [lCj 



•si = G11 + (d- 1)G 66 > 0, s 2 
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G11 — G12 — 2G66 , 1( ,n 
s 3 = • (19) 
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In terms of Every's parameters the inequalities (|17[) and 
HH1> read 



si > 0, s 2 < 1, 



(20) 



2D : (1 - s 2 + s 3 ) > 0, (s 2 - s 3 /2) > 0, 

(21) 

3D : (1 + 2s 2 ) > |4s 2 - 3s 3 - 1| , (10s 2 - 6s 3 ) > 1. 

(22) 

For both quadratic and cubic materials inequalities 
(|20| -([22 |) define stability regions in a form semi-infinite 
prisms (si > 0) with two distinct triangles in the base. 
Both stability triangles (STs) are l ying in the (s 2 , sn)- 
plane (Figs. QJand ffl Acoustic ([l(f, [TD] and [l2|) 
and elastic properties , Q of a particular initially un- 
strained cubic and quadratic material are characterized 
by si, which has dimensionality of stiffnesses, and two 
dimensionless parameters s 2 and S3 belonging to the ap- 
propriate stability triangles. Isotropic materials are rep- 
resented by points of the interval of the line S3 = be- 
longing to the suitable ST. For acoustic and elastic char- 
acteristics the parameter si is a scaling parameter. Re- 
gions of stability for symmetry systems other than cubic 
and isotropic as well as for the oblique symmetry sys- 
tem are defined in spaces of dimensions higher than 3. 
In the case of rectangular symmetry system the space of 
parameters can be effectively reduced to a 3D space [l]| . 



IV. PROCEDURE OF AVERAGING OVER 
DIRECTION OF CRYSTALLITES 

Following the assumptions made at the end of Sect. U 
to find the mechanical characteristics of polycrystalline 
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FIG. 1: Division of the stability triangle for monocrystalline 
quadratic elastic materials into areas of various auxetic prop- 
erties. 
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FIG. 2: Division of the stability triangle for monocrystalline 
cubic elastic materials into areas of various auxetic properties. 



materials one should calculate mean values of expressions 
defining them. Consider quadratic materials and a func- 
tion F of angle ip, then in agreement with definitions Q 



In the case of cubic materials one should consider func- 
tion F of angles a, ip and 9 introduced in Sect. Ill Al The 
mean value (F) of F(a, ip, 0) is defined as 



+7T pTT flit 

da I dip simp / d6F(a,(p,8). (24) 

- 7T JO JO 



Since we consider the model of randomly oriented 
crystallites, the averaged mechanical characteristics are 
isotropic. This means that if {v) > 0, one deals with an 
isotropic non-auxetics, whereas if (v) < 0, a polycrystal 
is a complete isotropic auxetic. We conclude that in the 
case of considered models of polycrystals the division of 
the stability triangles simplifies. 

After averaging v q (ip) = (i/ q (ip)/E q (ip)) E q (<p) we ob- 
tain 



(S2, S3) 



f (s2-l)(s 2 -S3~l) 

(S 2 + 1) (>2 - S 3 + 1) ' 



(25) 



V. DIVISION LINES FOR CRYSTALLINE AND 
POLYCRYSTALLINE QUADRATIC AND CUBIC 
MATERIALS 

Figs. 1 and 2 present the lines dividing the stabil- 
ity triangles into areas of various auxetic properties for 
monocrystalline quadratic and cubic materials. Since 
each part of STs has a finite area, the property of aux- 
eticity is quite widespread. 



Consider quadratic materials. The condition 
yields an equation of dividing line 

10s 2 + 3 



» 



3 S 2 



3s 2 



= 



(26) 



Since E q (n) ^ one can consider a different line de- 
fined by the condition (v q /E q ) = 



O2) = 3s 2 - 2 + yj3sl-4s 2 + 



dipF(<p). 



(23) 



(27) 



Lines ((26)) and (|27|) intersect at points s 2 = ±1 and on 
the line S3 = 0. Inspecting Fig. [3] we note that for both 
considered polycrystalline materials the area of complete 
polycrystalline auxetics is larger than for monocrystalline 
auxetics. This means that there exist cubic and quadratic 
auxetic materials that become complete auxetics in the 
polycrystalline phase. 

Graphene has hexagonal symmetry, i.e. it is an 
isotropic material. Its elastic properties depend on two 
stiffnesses, namely Cn and C%2- They were calculated by 
Falkovsky [3, Michel and Verberck fijl, and Kliten- 
berg et al [17| and measured (cf. [14[ )■ Points on the 
S3 = isotropy line standing for these data are depicted 
in Fig. [U It is seen that in both monocrystalline and 
polycrystalline 18] phases graphene is a non-auxetic elas- 
tic material. A point representing single ZnO monolayer 
with graphene-like structure [l9[ is also set forth. Two 
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FIG. 3: Comparison of division lines of ST for monocrys- 
talline and poly crystalline quadratic elastic materials. Points 
of area with lighter shade represent complete monocrystaline 
auxetics, whereas in the area with darker shade lie complete 
polycrystaline auxetics. 



FIG. 5: Comparison of division lines of the ST for monocrys- 
talline and polycrystalline cubic elastic materials. Points of 
area with lighter shade represent complete monocrystaline 
auxetics, whereas in the area with darker shade lie complete 
polycrystaline auxetics. 
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FIG. 4: Location of graphene on the isotropy line. The loca- 
tion of graphene-like ZnO is also indicated. Graphene refer- 
ences: 1-3 - [lfl, 4 and graphen exp. - Reference for 
ZnO: QJ|. Values of si (in 10 GPa) for graphene: = 56.2, 



= 58, 



(26) 



= 71, s 



(3) 



for ZnO: s 



(i) 



97.5, 



J 2 ) 



= 53.7, 4 4) 



120, s 



(exp.) 



= 150, 



65.1. Results of Klintenberg et al 



[171 ] for graphene are located between points 3 and 1. 



points representing single ZnO monolayer with graphene- 
like structure is also set forth. 

In the case of cubic materials numerical calculation 
gives the line of vanishing (v c ) = which is presented in 

Fig. m 



VI. MEAN VALUES OF MECHANICAL 
CHARACTERISTICS OF POLYCRYSTALLINE 
MEDIA 

In our previous papers we presented maps of E, G 
and v for quadratic materials Qand cubic p[. Since 
for monocrystalline media the mentioned characteristics 
are anisotropic, these maps refer to n and m directed 
along high symmetry crystalline directions. In the case of 
polycrystalline materials mechanical characteristics are 
isotropic. Their maps are presented in Figs. [S]and[7J 

Calculating (E q ), (G q ) and {v q ) we obtained 



(E g ) = 



4CJC L C M 



\J[2clcm + cj (cl + c M )] 2 - Cj (cl - cm)" 



= 1 - 



(Gq) = V^7/2, 



Ac L c M 



(28) 
(29) 



\J [2c L c M + cj (c L + c M )] 2 - c] (c L - c M ) 2 

(30) 

Table U summarizes the properties of mean val- 
ues of the mechanical characteristics and eigenvalues 
cj (I = J, L, M) on sides of ST. These properties align 
with calculated maps of (E q ), (G q ) and (v q ) (Figs. [5]and 
El). 



■5 




FIG. 6: Maps of (E q ), (G q ) and (v q ) for quadratic materials. 
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FIG. 7: Maps of (E c ), (G c ) and (v c ) for cubic materials. 



Generally, for various points of (S2, S3)-space, i.e. for 
various quadratic or cubic materials, the calculated mean 
values arc different. Since on the line S3 = both 
quadratic and cubic materials are isotropic, the maps for 
monocrystalline and polycrystalline materials should co- 
incide along this line. This requirement is fulfilled with 
a good accuracy. 



VII. CONCLUSIONS 



This paper has presented properties of the model of 
initially unstrained 3- and 2-dimensional cubic polycrys- 
talline solids. We have assumed that grains have the 
same volume and shape, and are randomly oriented. We 
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TABLE I: Values of (E q ), {G q ) and (v q ) on the sides of ST. 
J - top side, L - vertical side, M - bottom side. 



Side 


vanishing eigenvalue 


(E q ) 


(G q ) 




J 


C.T 







-1 


L 


CL 








1 


M 


cm 








1 



have shown that some monocrystalline auxetics become 
in the polycrystalline phase complete auxetics. Gener- 
ally, the averaged mechanical characteristics follow from 
the weighted integration over all orientations. Weighting 
function is the orientation distribution function which de- 
scribes the texture of the polycrystal considered (eg. fibre 
textures). These more general models of polycrystalline 
solids certainly deserve further studies. 
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Erratum 

We noted three mistakes in article by Jasiukiewicz et 
al, physica status solidi vol. 245, issue 3, pp. 562-569: p. 
564 formula for tq should read: 



+ 2Si2 — 5*11 — S[ 



S 



16 



S 2 



r G 



The first of Eqs. (13) should read: 

E^ 1 = Sn cos 4 ip + S22 sin 4 ip + 2Siq cos 3 ip sin ip 

+2S26 cos (p sin 3 ip + (2512 + 566) cos2 ^ sm2 V' 
p. 565 the definition of K should read: 

K = (Cn - C12 - 2C 66 ) 12. 



